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Abstract. The classical version of the three wave interaction models the creation 
and destruction of waves; the quantized version models the creation and destruction 
of particles. The quantum three wave interaction is described and the Bethe Ansatz 
for the eigenfunctions is given in closed form. The Bethe equations are derived in 
a rigorous fashion and are shown to have a thermodynamic limit. The Dirac sea of 
' negative energy states is obtained as the infinite density limit. Finite particle/hole 

, excitations are determined and the asymptotic relation of energy and momentum is 

' obtained. The Yang- Yang functional for the relative free energy of finite density exci- 

tations is constructed and is shown to be convex and bounded below. The equations 
of thermal equilibrium are obtained. 
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1. INTRODUCTION 

The creation and destruction of particles is a fundamental process in high energy 
physics. Among other things, the presence of such processes implies that the en- 
tries of the scattering matrix between states with different numbers of elementary 
excitations (or physical particles) can be non-zero. 

Integrable systems provide solvable paradigms for various physical phenomena: 
e.g. for lattice statistical models, the algebraic decay of correlation functions at 
a critical point may be observed. The creation and destruction of particle pairs 
is modeled by the quantum three wave interaction (3WI), which is obtained by 
quantizing a particular reduction of the = 3 case of the AKNS hierarchy as 
in Wadati and Ohkuma^, or by taking the continuum limit of the GL{3) lattice 
magnet as in Kulish and Reshetikhin^ . The result is an interacting many body 
system which does not conserve the number of particles: the interaction permits 
two distinct types of particles to combine to form a third, and also the reverse 
process. 

The quantum 3WI has been considered in three forms. In one, all particles are 
bosons, in the "Fermion I" model two fermions interact to form a boson, and in 
the "Fermion II" model a boson and a fermion interact to form a fermion. In this 
paper we consider the Fermion I model. We find the associated first quantization 
problem, obtain the eigenfunctions of the Hamiltonian in closed form, and derive 
the Bethe equations rigorously from the requirement that the wave functions be 
periodic. 
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Because the model is first order in the derivatives, the energy is not bounded 
below. Nevertheless, as in Dirac's model of the electron, a stable ground state may- 
be constructed by filling all negative energy states. One can then consider all states 
with reference to the filled Dirac sea of negative energy states. This approach is used 
in Korepin"^ and in Bergknoff and Thacker^ to calculate the excitation spectrum of 
the massive Thirring model from the Bethe ansatz states. 

We prove the existence of a thermodynamic limit as the length of the interval 
tends to infinity with constant density. The Dirac sea is obtained in the infinite 
density limit. Finite particle/hole excitations from the Dirac sea correspond to 
the solutions of certain integral equations. Such excitations have positive relative 
energy, and we determine the asymptotic relation between energy and momentum. 
Following the treatment of a one-dimensional boson gas with delta-function interac- 
tions by Yang and Yang^, we derive a functional for the relative free energy of finite 
(local) density perturbations from the Dirac sea. We prove that the functional is 
convex and bounded below and obtain the equation for thermal equilibrium. 

Kulish and Reshetikhin^ used the quantum inverse scattering method (also called 
the algebraic Bethe Ansatz) to investigate GL{N) invariant transfer matrices on 
a one dimensional lattice. They obtained quantized N wave models as a formal 
continuum limit and found explicitly the Bethe equations for the bosonic 3WI. 
Wadati and Ohkuma^ also considered the bosonic case. Ohkuma and Wadati^ 
treated the Fermion I and Fermion II cases and gave an algorithm for constructing 
the Bethe wave functions. 

Ohkuma^ investigated the Fermion I quantization of the three wave interaction, 
wrote down the Bethe equations, and calculated the free energy and excitations 
from a ground state. However a momentum cut-off is introduced in Ohkuka^. The 
results depend on the cut-off and do not have a finite limit when the cut-off is 
taken to infinity. By constructing all physical states with reference to a Dirac sea 
we obtain that are finite and are independent of any cut-off. The methods provide 
a basis for calculating the free energy of other integrable models whose spectrum 
is unbounded below. 

The quantized field of the 3WI may be rewritten as a first quantization for 
a many body problem. Such a reduction was carried out by Lieb and Liniger^ 
for the quantization of a Bose gas - essentially the quantization of the nonlinear 
Schrodingcr equation. They obtained a Helmholtz equation in the interior of the 
fundamental domain, together with boundary conditions at the interfaces associated 
with the interactions; see also Korepin, Bogoliubov, and Izergin^. A novel feature 
of the 3WI is that the system of equations links functions with different numbers 
of variables, corresponding to the fact that individual particle numbers are not 
conserved. 

As in Lieb and Liniger^ we solve the differential equations in certain fundamental 
cones C'^''^'P and then extend the densities to all space, first by anti-symmetrizing 
across the boundaries of C'^''^'P and then by periodicity. The differential equations 
for the densities force discontinuities where fermions of different types interact. 
However like fermions do not interact, so to avoid discontinuities on the boundaries 
of the fundamental cones we impose a regularity requirement, (2.5). This require- 
ment forces one to anti-symmetrize the basic solutions of the differential equations 
over momenta. This in turn is fundamental to obtaining the Bethe equations, as 
we show by example in §5. 

T„ ceo A „i ;„ „ J- — J„j-„;i ^-^ ;„ ] ] 
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Ohkuma^, a derivation of an algorithm for the eigenfunctions of the 3WI, and go 
on to obtain the explicit form of the eigenfunctions (cf. Theorem 5.1). This form 
generalizes the Bethe Ansatz (Bethe^°, Lieb and Liniger^, Lieb^^, Yang and Yang^^) 
to systems in which the eigenvectors are superpositions of states having different 
numbers of particles. In §5 we derive the Bethe equations from the requirement of 
periodicity. 

In §6 we reformulate the finite Bethe equations as a system of integral equations 
and prove existence of a certain thermodynamic Umit at finite density. The Dirac 
sea is then obtained as an infinite density Hmit. In §7 integral equations for the 
particle/hole excited states are obtained, and energy /momentum asymptotics are 
found. Finite density perturbations of the Dirac sea are constructed in §8. In §9 
we investigate the (relative) free energy functional. The formal Euler-Lagrange 
equations for its minimum are the same as those proposed by Ohkuma^ if one 
imposes a cut-off. Without a cut-off they have no strict solutions, but there is 
a natural interpretation which does yield meaningful equations for the thermal 
equilibrium state. 

2. THE THREE WAVE INTERACTION 

The Hamiltonian of the classical three-wave interaction is 



where '^j are complex valued functions and f denotes complex conjugation. The 
constants Cj are wave speeds, and 7 is a coupling constant. Zakharov and Manakov-*^^ 
found a Lax pair for the classical system and Beals and Sattinger-^^ gave a set of 
action-angle variables. However, intcgrability of the classical system will play no 
role in our analysis of the quantized model. 

The system is formally quantized by taking the , \1/ j to be creation and anni- 
hilation operators on a Fock space, f denoting Hermitian adjoints. We take 'ifo{x, t) 
to be a boson field operator, and ^^i, to be fermion fields. Thus the operators 
are required to satisfy the equal time commutation ([ , ]-) or anti-commutation 
([ , ]+) relations 



The Hermitian operator T represents the kinetic energy, with distinct wavespeeds 
Cj for each field. The interaction term formally anihilates the two particles with 
wavespeeds ci and C2 (particles of type 1 and 2) and creates one with wavespeed 
Co, while V destroys a type particle and creates types 1 and 2. 

We assume throughout that the speeds of the fermions bracket the boson speed: 
(ci — co)(c2 — Co) < 0. In addition, we make the simplifying assumption that 




■^j 7(^'T^i^'2 + *^*l*o) dx 



(2.1) 



=T + 7(W + y), 



{x,t), {x\ t)]+ =0, [-^j (x, t) , {x\ t)]+ = 5,,,S{x - x') , 
[^o{x, t), ^o{x', t)]- =0, [^o{x, t), ^l{x', t)]_ = 6{x - x'), 
[*o(a;,t),*,(a;',t)]_ =0, [*J(a;, t), *,(a;', t)]_ = 0, j,/ = 1,2. 



(2.2) 
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where a is an arbitrary constant. Without loss of generahty we assume throughout 
this paper that a > 

The Hamiltonian is defined on a Hilbert space Ji consisting of the tensor product 
of a boson and two fermion Fock spaces. There exists a non-zero vector |0) in the 
Fock space, called the pseudovacuum, such that 

*,(a;,t)|0) = 0, j = 0,1,2. 

Define particle number operators 

Nj = J dx^]^j, j = 0,1,2. 

The coordinates of particles of types 0,1, and 2 are denoted by Xj, yj, Zj respec- 
tively. Let 

^l{x) = ^l{xi,...,Xm) = ^l{xi) ■ ■ ■ ^l{Xm) 

and define ^\{y) and ^2i^) similarly. Then set 

*t(x,y,^) = vE't^x)\E'l(2/)\E'^(^), xeR'^^yeR^, z e RK 
Fundamental states in Ti, are given formally by 

^\x,y,z)\Q), 

which are eigenstates of Nq,Ni, and N2 with eigenvalues m,n, and p respectively. 
A general state vector in this subspace may be written in the form 

|'0m,n,p) — j j T / fmnp 

{x,y,z)'l>\x,y,z)\Qi)dxdydz, (2.4) 

Ul'.n'.p'. J J^m,n,p 

where fmnp belongs to C2{R^'"''^) and dx = dxi dx2 ■ ■ ■ dxm, etc. Since ^'^(x, y, z) 
is symmetric in the Xj and anti-symmetric in the y^ and zi, the density function 
fmnp may be assumed to possess these same symmetries. 

If a density fmnp is associated with an eigenstate, then, as shown in §4, it must 
be discontinuous across the hyperplanes y^ = zi, which we call the interaction 
manifolds. It is natural to assume that these are the only discontinuities. We 
therefore assume continuity across the hyperplanes yk = yk+i, provided yk ^ zi 
for all zi, and similarly with respect to the z coordinates. Continuity plus anti- 
symmetry across these hyperplanes implies: 

Regularity condition : fmnp vanishes as yk yk±i if 

yk±i 7^ zi for all I and similarly as zi zi±i 7^ yk, all k. (2.5) 

By symmetry, we may replace the integration in (2.4) by integration over the cone 

(lrn,n,p g^yg^ by 

X\ <X2<--- < Xm] yi < y2 < ■ ■ ■ < yn, Zi < Z2 < ■ ■ ■ < Zp, 

and write 

l^m.n.v) = I fmnv{x,y,z)'^Hx,y,z)\Q)dxdydz. (2.6) 
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Prom now on, we write the state vectors in the form (2.6) and assume that the 
densities fmnp are supported on the cone C'^''^'P and satisfy the regularity condition 
(2.5), which implies that the antisymmetric extensions of the fmnp are continuous 
across the collision manifolds yj = yk and Zj = Zk of like fermions. 

The classical version of the Hamiltonian (2.1) is integrable and possesses an 
infinite number of conservation laws: Zakharov and Manakov^^. Three important 
conserved quantities for each version are 



N = J 2*J*o + + *2*2 dx 



i=o 

In the quantized version operators N and D give conserved particle numbers, and 
P is the momentum operator. Indeed it is not difficult to show from (2.2) that 
N,D,P commute with H. 

3. REDUCTION TO FIRST QUANTIZATION 

The quantum field problem can be rewritten as a first quantization problem for 
a many body problem. 

Denote the subspace of states \iprn,n,p) by Tmnp- These are eigenstates of N 
and D with eigenvalues N — 2m + n + p and D = n — p respectively, but are not 
eigenstates of the Hamiltonian H, since and V create and destroy particles. 
However, since N and D commute with H, the quantum numbers N and D label 
a subspace that is invariant under and D. We denote this space by J^n,d- 

We now fix N and D. Within Tn,Di densities can be labeled by m alone and 
we denote the corresponding subspace by Tm- Note that both of ± D must be 
non-negative even integers. When there are no bosons present, (m = 0), there 
are {N ± D)/2 fermions of types 1 and 2 respectively. The maximum number 
of bosons occurs when the maximal number of fermions have combined, yielding 
M = min {N ±D)/2= |(7V - \D\). 

The operators T,V\ and V act on the subspaces J^m as follows: 

T '■ J'm ^ mi '• -F m ^ F m+li ^ '■ Fm ^ F m — 1- 

An eigenstate of the Hamiltonian in J-n,d is a sum of vectors iV'm) £ Fm-, and the 
eigenvalue equation becomes the coupled system 

{E-T)\^Po) =7^|Vi) 

{E - T)\'^rn) =lV^\'^m-l) + lV\'^m+l) (3.1) 

/ IT' rr^W^I. \ T/tl„/. \ 
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since V|V^o) = V^i^M) = 0. 

Let fm be the density corresponding to the vector iV'm)- Equations (3.1) may be 
converted to a system of partial differential equations by calculating the action of T, 
and V on the densities fm- The calculation of the action of T is straightforward, 
while the calculation of V and V"^ is somewhat more involved; details are given in 
Beals, Sattinger, and Williams^^. The action of T is: 

T\iPm)=T [ fm{x,y,z)-^\x,y,z)\0) dxdydz 

J Qm,n,p 

{Xmfm{x,y, z))'^'^ {x,y, z)\0) dxdydz 



where 



j=o fe=o ^y'' 

We use to denote the action on densities as well as the action on vectors. Then 

m+l,n,p 

v^fm= Yl {-ir^'^^+''+'fmi£j;{y,^j,k);{z,xj,i)), (3.5) 

j,k,l=l 

where {x,y,z) e cm.+i,n-i,p-i ^^^^ 

— ("^1) • • • ) Xj—\^ a^j-i-i, . . . x^f^ G C , 
(y, Xj,k) = {yi,..., yk-i,Xj, yk,---, yn-i) e C", etc. 

Remarks. 1. Since fm is supported on C'^'^'P^ fm{xj', {y, xj, k); {z, Xj,l)) = unless 
Xj e {yk-i,yk) n {zi-i,zi). Hence 

{V^fm){x; y;z)^Y. (-i)"^"'^+'^+'V-(^,; {y, Xj, k,); {z, xj, 1^)), (3.6) 

i=i 

where kj and Ij are uniquely determined by the constraint 

e ivkj-i, ykj) n {zi.-i, zi.). (3.7) 

2. In (3.6) we interpret the value of fm on the interaction manifold Xj = yk^ = zi. 
as the average of its values as yk^ —>■ zi.±. 

3. The action of V'^ can be thought of as gluing a type 1 fermion at and 
a type 2 fermion at zi to form a boson at Xj. The coordinates are relabeled to 
take into account the shift in numbers of variables, summed over all possible pairs. 
Similarly, the action of V can be thought of as the decomposition of a boson into a 
pair of fermions, summed over the bosons. In fact the operation of V on densities 
is given by 

n+l,p+l 

{Vfm){x,y,z)= i-'^r^'^^^''^'fm{{x,ykJk);yk;zi)5{yk-zi). (3.9) 
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where jk is uniquely determined by the constraint i/k G {xj^^-i, xj^). 

We illustrate in some low order cases. For 2m + n+ p — N — 2, fo — 

fo{y,z), fi = fi{x): 

{V^fo){x) = -fo{x,x), Vfo = 

^Vi=0, iVh)iy,z) = hiy)5iy-z) 

For 2m + n+p = iV = 4, fo = fo{yi,y2, zi, Z2), fi = fi{x,y, z), f2 = f2{xi,X2): 

{V^fo){x, y, z) =fo{x, y, x, z) - fo{x, y, z, x) - fo{y, x, x, z) + /o(2/, x, z, x), 

{V^^ fl){xi,X2) = - fl{xuX2,X2) - fl{x2,Xi,Xi), 

VU2 =0; 
V^/o=0 

{Vfi){yi,y2,zi,z2) =fi{yi,y2,zi)5{yi - Z2) - fi{yi,y2,z2)S{yi - zi), 

- fiiy2,yi,zi)S(y2 - Z2) + fiiy2,yi,Z2)S{y2 - ^i), 
{Vf2){x, y, z) =[f2{y, x) + ^(a;, y)]6{y - z). 

4. THE BETHE ANSATZ 

The differential equations for the densities, supplemented with the regularity 
condition (2.5) at the boundaries of the cones C"*'"'^, describe the spectral problem 
in first quantization. 

From the results of §3 the differential equations for the densities are 

Xofo + Wfi^Efo, 

Xmfm + lVfm+1 + iV^ fm-1 = Efm (4.1) 

Wadati and Ohkuma-^ gave an algorithm for constructing the general eigenstates. 
In the next two sections we find these eigenstates in closed form, in two steps. In 
this section we construct "local solutions" of equations (4.1); on each of the cones 
Qmnp ^ multiple of an exponential solution of the associated differential equation. 

Because the operator V introduces delta function singularities on the interaction 
manifolds, solutions must have discontinuities across these manifolds. However the 
local solutions violate (2.5) and thus they have additional discontuities if extended 
by anti-symmetry. We denote the densities of the local solutions by hm to distin- 
guish them from the regular densities, which will satisfy (2.5). We obtain these 
regular densities in §5 by taking linear combinations of the hm- anti-symmetrizing 
over the momenta. 
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Two low order cases illustrate the main features of these equations. For N = 
2, D = 0, equations (4.1) are 

*ci^ + ^C2^ + ihi{y)5{y - z) =Eho{y, z) (4.2a) 
dhi 

ico-^ - ^ho{x,x) =Ehi{x) (4.2b) 

where ho — ho{y, z), hi — hi{x). In the first equation the interaction term vanishes 
off the interaction manifold y = 2; so we take the solutions to be multiples of an 
exponential in each of the open sectors y < z and y > z oi 



ho{y, z) = 



ex.pi{r]y - Cz), y < z, 
ex.pi{r]y - (z), y > z, 



where E = —{ciT] — C2C)j and (pi and 02 ^-re constants to be determined. 

In (4.2b) we take ho{x,x) on the interaction manifold to be the average of the 
values of ho{y, z) from the left and right. Hence the second equation is 

(,<,A-i5).,W^.^e'<-««. 

We seek a solution of the form h\{x) = ^e*^^~^^^, so that 

hi{x)=g{v,C)^4^e^pi{v-0x, 9{v,C)- ^ 



(ci - co)r] - (c2 - co)C 
By (2.3), 

5'(^,C) = ^^-:t:=^(^ + C)- 

a{ri + C) 

We use this form of g from now on. 

Finally, we determine the ratio 4>i/(j)2 by integrating (4.2a) across the interaction 
manifold y = z. We change variables, letting d = y — z, s = y + z, and obtain 

«(ci - C2)^j + i{ci + C2)^^ + lhi5{d) = Eho, 

, ,'s + d s - d\ , , f s + d 

ho =ho —z—, —z— , hi = hi 



2 ' 2 

Integrating across the jump d = {y = z = x) we get 

i{ci -C2)(02 -</'i)expia;(?7-C) +7^(^ + ^^-^expia;(77-C) = 0; 

hence 

^ ^ 2z(ci - C2) - 'jg 
4>\ 2z(ci - C2) +75'' 

Using our assumption (2.3) we may write this in the simplified form 
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where 

UJ = 



Setting 05 = 1, we have 



ho{y,z) 



expi{r]y - (z), y < z; 

e{r] - C) ex.pi{r]y - (z), y > z; 



and 



hi{x) = i^(?7 + C) (1 + 0{r] + 0) expix{r] - ()■ 



For N = 3, D = 1 equations (4.1) are 

. dho dho dho 

ici ^ h ici h iC2 EhQ [yi ,y2;z) 

oyi dy2 dz 



=7(^1 (yi; 2/2)5(2/1 - z)- hi{y2]yi)d{y2 - z)) (4.3a) 

(4.3b) 



^co^ + ici^ - Ehi{x;y) = -f{ho{y, x; x) - ho{x,y;x)), 



with ho = ho{yi, j/2; z) and hi = hi{x; y). Arguing as above, we obtain 

{expz(?7iyi+?72?/2-C^), yi<y2<z; 

0{V2 + C) expi{T]iyi + ?72?/2 - C^), yi < z < 2/2; 

9{r]i + C)9{r]2 + C) expi{r]iyi + 7722/2 - C^), z < yi < y2] 

and 

1 r (1 + e{r]2 + C)) 9{t12 + exp i{{r]2 - 0^:^ + VW) y < x; 



hi{x; y) = 

2 I -0{ri2 + 0(1 + 9{Vi + OgiVi + C) exp i{{rii - Qx + 7722/) x <y. 
The General Case 

The two preceding cases show that when y^ and zi are transposed, (with y^ 
moving to the right of zi), Hq is multiphed by 9{'r]k + Ci)^ where rjk, —Q are the 
momenta of the corresponding fermions. Since the differential equations are local 
and involve only pairwise interactions, this situation extends to the general case. 

We begin by introducing some notation and conventions. For fixed N and D, 
(^rn,n,p uniquely determined by m, so from now on we denote it by C^. The 
domain is a union of 

{m + n + p)! 
m\n\p\ 

open sectors determined by the interspacing of the particles of type 0,1,2. We label 
them by m + n + p-tuples cr = (ei, €2, . . . , e^+n+p) where Cj — 0, —1, or +1 if the 
particle in the j^^ position along the real line is of type 0,1, or 2 respectively. For 
example 

a = (1, 0, -1, 1) ^ {zi <xi <yi < Z2}. 
The boundary of C"^ consists of the interaction manifolds yj = Zk together with 
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A transposition consists of interchanging a -1 with a +1 to its immediate right, 
e.g. (1,0,-1,1) I— > (1,0,1,-1). Each transposition corresponds to interchanging 
an adjacent pair of coordinates yk, zi; thus to every interaction manifold there is a 
corresponding transposition. The sectors of C"^ are partially ordered by the number 
of transpositions; thus, for example (1,0,-1,1) < (1,0,1,-1). By definition, no 
transpositions occur between -I's and I's separated by one or more zeroes. 

Since Vhi is supported on the interaction manifold, ho = ho{y; z) satisfies 

fri dyj ^ dzk 

in each open sector o" C C°. 
Let 

cro = {yi < ■ ■ ■ < yn < zi < ■ ■ ■ < Zp} (-1, . . . , -1, 1, . . . , 1). 

We take 

(n p \ 

^ VjVj - ^ CkZk , (y; z) e ao 
3=1 k=i J 



with 



E = -ci5^?7,+C2 5^Cfe- (4.5) 

j=i k=i 



The solution in a general sector a G C° is 



ho = 0° expz r)jyj - CkZk , y, z e a, (4.6) 
\i=i fe=i / 

where the constant (/)° is uniquely determined by cr, and is given in (4.8) below. 

Lemma 4.1. Let a < a' differ by a single transposition of fermions of types 1 and 
2 at yj, Zk with momenta rjj and — Cfe respectively. Then 

0°, = ^(r/,+a)0°, (4.7) 



Remark. For a' < a the phase factor is 

nri,+Ck)]-'^e{-r^^-Ck). 

Thus within a fixed particle number sector these conditions on the wave function 
are of the usual Bethe Ansatz type; the ratio of phases that differ by a permutation 
of neighboring particles is a function of the momenta. 

Proof. All the interactions are pairwise and local, so the same computations for the 

case N = 2, D = extend to the general case. The ratio of the constants (p'^, (p^, is 
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Each cr e C° is obtained from ctq = (—1, ■ ■ ■ ■, —1, 1, • • • , 1) by a unique set of 
transpositions {yj,Zk), which we denote by T^. Thus 



n (4.8) 



Now consider the higher order components of the density: h^n, m > 1. The 
sectors of are characterized by {m + n + p)-tuples a with m zeroes. We assign 
momenta to a as foUows. Regard each in u as a pair (—1,1). Starting from 
the left, assign the -I's momenta 771,772, • • • ; the +l's momenta — Ci, — C2, ■ • ■ ; and 
assign the O's the sums of momenta rjj — C/-, with the ry's and C's taken in their 
normal order. For example, 

(-1,0,1,0) ^ (r/i,r/2 - Ci,-C2,^3 - Cs); 
(-1,1,1,0,1,0,-1)^ (r7i,-Ci,-C2,772 -C3,-C4,^3 -(5,^4). 

For fixed locations of the zeroes in a the basic exponential solution is defined by 

(m n p \ 

^riVir - Cjr) + ~ Y ^^^^^ I (^•^) 

r=l s=l t=\ ) 

where 

V (ir) = # 3;'s and j/'s (2;'s) to the left of Xj.-, inclusive of x^'-, (4-10) 
O'S ih) = # a;'s and j/'s (2;'s) to the left of ys, (zt), inclusive of j/s, (zt). 

Remark. From (3.7), ir — r + kr — 1, jr — r + Ir — 1- 
For example 

e^(a;; y; 2;) =expza;i (0^1(772 - Ci) + 2^2(773 - C3) + ViVi - ^1(2) , 
a= (-1,0,1,0); 

e^(a;; y; z) = exp z (a;i (7/2 - C3) + 3^2(773 - Cs) + ^1^71 + 2/2774 - ^iCi - ^2(2 - ^3(4) 
a =(-1,1, 1,0, 1,0,-1) 

We now describe the general case. For any a e define ^o- by 

m 

9CT=YI gUr-, 9ir3r = 9{Vir + 0. ) , (4-12) 

r=l 

where v, jr are given in (4.10). For example, for (—1, 0, 1, 0), ii = 2, ji = 1, Z2 = 
3, j2 = 2, and 

9a = 921932- 

The components of are determined by the location of the zeroes and the 
numbers and orders of +l's and —I's between successive zeroes. The factor g^r is 
constant on each connected component of C^. For a e C"^, we define 

C = f[(-l)nin^)+kr+lrl±^in^r±(iA^ (4.13) 

where S G C'' is obtained from a by replacing each in cr by (—1, 1), and (f)j. is 
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Theorem 4.2. The basic solutions of equations (4-1) are given by 

hm{x,y,z) = ip:^e'^{x,y,z), (4.14) 

where 




and the 0^ are given by (4.8), (4-13); the g^ by (4-12); and the by (4-9). 

The proof is given in Beals, Sattinger, and Williams^^. In addition to Lemma 
4.1, it relies on the identity 

(_l)n(m)+fc.+/. Ifr 1^ = 0^+1^ 1 < r < m + 1 

where cr^ and a'^ are obtained from a by replacing the r-th zero by (—1,1) and 
by (1,-1) respectively. This condition insures compatibility between the phases in 
different particle number sectors; it has no analogue in the usual Bethe Ansatz. 



5. THE BETHE EQUATIONS 

In §4 we constructed basic solutions km of equations (4.1). As noted in §2 we 
want to impose the regularity condition (2.5) that eigendensities have no unneces- 
sary discontinuities, i.e. that their anti-symmetric extensions be continuous across 
the collision manifold of like fermions. To obtain this property, we antisymmetrize 
the basic solutions over the momenta. 

Permuting the momenta of the basic solution of (4.4), (4.5) gives another solu- 
tion. The general solution of (4.4, 4.5) is a linear combination 

(n p 
j=i k=i 

where J, K are permutations of 1, . . . , n and 1, . . . ,p respectively. In particular, 

/o(y; z) = ^sgn(J) sgn(i^) eo(J, K), 

J,K 

satisfies Xq/q = Efo in 

(Jo ^ {yi < ■ ■ ■ < Vn < zi < ■ ■ ■ < Zp}, 

and satisfies the regularity condition (2.5). 
The density in the other sectors is given by 

My;z) = J]sgn(J)sgn(i^)0°(J,i^)eo(J,i^), (5.1) 

J,K 

where {y;z) e a, and 
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The corresponding densities for higher m are consequently 

U{x; y;z)=J2 sgn( J) sgn(K) C( J, K) g^{J, K) e^{J, K), (5.2) 

J,K 

{x;y;z) e a, 

where e^(J,K), ga{J,K) and (j)^{J,K) are obtained from (4.9), (4.12), and (4.13) 
by making the substitutions j, k ^ Jj, Kk- For example, to obtain densities which 
vanish when yi = y2 z m. the case N = 3, D = 1, we antisymmetrize the base 
function over the momenta: setting 

eo(j, k) = expi{r]jyi + r]ky2 - C^), 

we take 

r eo(l,2) -eo(2,l), yi < ?/2 < z; 

/o(yi, 2/2, z) = I e{7]2 + eo(l, 2) - 0(?7i + eo(2, 1), ?/i < 2; < j/a; 

[ ^(771 + 0^(772 + C)[eo(l, 2) -eo(2,l)], 2 < yi < ya- 

Now fo{yi,y2, z) vanishes when j/i = ?/2 7^ -2; but we can not assert anything for 
the case yi < z < y2, since yi and ?/2 cannot coalesce without equaling z. The 
corresponding density /i is also obtained as a superposition: setting 

k) = g{r]k + C) (1 + 0{r]k + C)) ewKVk - C)x + VjV) 

we find 

h{x, y) j ^ ^^^^2, 1) + e^rj^ + C)w{l, 2), x<y. 

Theorem 5.1. The antisymmetrized wave functions (5.1) and (5.2) satisfy the 
regularity condition (2.5) and thus determine a basis of eigenstates for the Hamil- 
tonian (2.1). 

Proof. Consider a pair (j/^, j/fe_|.i) (the proof is the same for contiguous 2;'s) and 
assume there are no intervening x' s or z' s. Group the terms in the sum (5.1) into 
pairs corresponding to J = (. . . , j, j + 1 . . . ) and J' = {. . + . . .). At yj = yj+i 
we have eo(J, K) — eo{J', K), and the coefficient of this common value is 

</.°(J,K)-</.°(J',K). 

For any given pair (j, /c) either yj < yj+i < Zk or Zk < yj < yj+i- In the first 
case, neither ( j, /c) nor (j + 1 , fc) belong to , so neither of the factors 9 {rjj. — ) , 
^{VJj+i ~CKk) belong to the products 0° (J, K), (f)^{J', K). In the second case, both 
(j, k) and {j + 1, k) belong to T^-, so both factors belong to both the products, and 
0^( J, K) = (j)^^{ J' ,K). Therefore the coefficient of eo( J, K) = eo( J', K) vanishes in 
either case, and the sum is zero. This establishes the result for /q. 

We proceed in the same way to verify that the higher order densities fm-, m>l 
vanish when yj and j/j+i coincide. In view of (4.16), we must show that 
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From (4.12) the constants g^{J,K) depend only on the momenta associated with 
the bosons (the zeroes of u), and so are unchanged under any permutation of the 
fermion momenta. The same applies to the product in (4.13): it depends only 
on the boson, not the fermion momenta. Therefore the issue for the higher order 
densities reduces to the proposition (f)j.{J, K) = (f)j.{J',K), which was established 
in the first step of this proof. □ 

Given the fm defined on the cells 

{ < a;i < • • • < a;^ < -C', < yi < • • • < y„ < L, < zi < ■ ■ ■ < Zp < L }, 
we extend them to the hypercube < Xi, yj, Zk < L by 



where Mx, Jy, Kz are the unique permutations of the variables j/, z for which 

XM^<■■■<XM^, yji<---<yj„, ZK^<---<ZK^. 

We now impose periodic boundary conditions on the wave functions hm- These 
periodicity conditions lead to constraints on the momenta, known as the Bethe 



hm{x\ y; z) = sgn{J)sga.{K)fm{Mx, Jy, Kz) 



(5.3) 



equations. 

For example, in the simplest case N = 2, D = (cf. §4) we obtained 




The periodicity conditions 



fo{0,z) 



fo{L,z), 



/o(2/,0) = /o(2/,L) 



lead immediately to the Bethe equations 



l = e{r] + Qe'"'^ 



For the case n — 2, p — 1, the conditions (5.3) lead to 



fo{yi,y2,0) = fo{yi,y2,L) 



fo{0,y2,z) 



fo{y2,L,z). 



From the definition of /o in this case we get, respectively. 



dim + C) dim + (exp iimvi + my2) - exp z(?7iy2 + V2yi)) 

= exp z(r/iyi + r]2y2 - LQ - exp iir]iy2 + mVi - ^C) 



and (say for y2 < z) 



exp i(?72y2 - C^) - exp i(?7i|y2 - C^) = 



-Oim + C) exp z(?7i2/2 + V2L - (z) 
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From these two equations we get the Bethe equations 

exp iLC = e{rii + C)0iV2 + Cy, 
1 = divi + C) exp ir]iL, 1 = 9(r]2 + C) exp ir]2L. 

All other periodicity conditions follow from these. 

Remark. The second example shows that one must anti-symmetrize in momenta 
to achieve periodicity. For example, suppose N — 3, D = 1, and < 1/2 < z. 
The condition ho{0,y2,z) — —ho{y2,L,z) for the basic wave function obtained in 
§4 becomes 

/io(0, 2/2, z) = exp i(?722/2 - C^) = -/io(2/2, L, z) = e{r]2 + C) exp i{r]iy2 + mL - (z), 
leading to the contradiction 

exp ir]2y2 = -0{r]2 + C) exp i{r]iy2 + ^L). 

Theorem 5.2. Let 

p 

e^^'=^=J]^(0 + %) k=l,...,n, (5.4) 

1=1 

n 

e<^L ^l[e{f^k + Ci) / = (5.5) 

k=l 

Then the densities km defined by (5.1), (5.2), and (5.3) are L-periodic in each 
variable. 

The equations (5.4), (5.5) are known as the Bethe equations. The proof is given 
in Beals, Sattinger, and Williams^^. 

6. CONSTRUCTION OF A DIRAC GROUND STATE 

A solution of the Bethe equations gives an eigenvector of the Hamiltonian with 
energy 

n p 

E = -ciY,Vk + C2j2Ci^ (6.1) 

k=l 1=1 

where n and p denote the maximum number of particles of types 1 and 2 for states 
in the space !Fn,Di i-e. n= ^{N + D), p = ^{N — D); cf. (4.5). The energy is not 
bounded below. We construct a stable ground state by filling all negative energy 
states. 

Taking logarithms in the Bethe equations (5.4), (5.5), we get 
2-K ^1 

Vk- j^Ik^'^Tj^^ogeirjk + Ci), k=l,...,n (6.2) 

1=1 

27V "1 
- = 5^T^log^(r7fe + 0) l = l,...,p 
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where 

6'(A) = —, CO = (6.3) 

and Ik are distinct integers, as are J/. 
Taking 

-log^(A) = 7r-2tan-^-, (6.4) 
i (jO 

we write the Bethe equations in the equivalent form 

2^ r ^2 'A _i r7z + ((,ri\ 
?7fc = — >^tan , (6.5) 

«=1 

n 



■ 2 j:tan-^^^+^' 

A:=l 



L a; 



where the 1^ are distinct half integers if p is odd, distinct integers if p is even, and 
similarly for the Ji in relation to n. We assume for convenience that both are even 
and, without loss of generality, that 

II < /2 < • • • < /n; Jl < J2 < • ■ • < Jp- (6.6) 

These conditions imply that there is a solution to (6.5) for which the rjk and the 
C,i are themselves strictly increasing. Indeed equations (6.7) are the variational 
equations for the functional 



\A;=1 1=1 J \k=l 1=1 J k,l=l 



2 



where 

.A 



12 s 
M(A) = - tan"^ - ds. 



For convenience in this section we take iv — 1. (This amounts to fixing a length 
scale, since cu has the dimension 1/length.) 

Theorem 6.1. The Bethe equations have at least one solution for each set of 
distinct integers Ik and distinct integers Ji. 



Proof. The term M(A) grows only linearly in its argument, so the quadratic terms 
dominate T, as r], ( ^ oo. It follows that S is bounded below and tends to infinity 
for large values of its argument. It therefore always has a minimum. □ 



Uniqueness of the solutions is a subtler question. For large values of the densities, 
S is not convex, so in principle it could have multiple critical points. In the limiting 
cases of interest, however, we will obtain uniqueness results. In fact we use the 
finite case to motivate the limiting cases. Therefore we save space and smooth 
the exposition by stating most of the finite case results without proof and leaving 
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We now turn to the construction of a stable ground state, or vacuum state. For 
the ground state we assume that D = (hence n = p = N/2) and we consider the 
case when the velocities of the fermions have opposite signs: 

C2 < < ci. (6.7) 

Thus ?7 > and C > mean negative energy for the corresponding particles. Our 
procedure will be to fill the available negative energy states from the top down by 
choosing consecutive integers Ik, Ji such that the solution {rjkXi} is positive but 
minimally so; we then take the infinite length and infinite density limits. The 
equations (6.5) are symmetric, so we take the data — Jk and look for a solution 
with rjk = Ck- Thus we want 

Ik = Jk = k - Kn (6.8) 
where is an integer to be chosen. 

Theorem 6.2. There is a largest integer Kn = 0{n) such that (6.5) with data 

(6.8) has a positive solution. This solution is unique and satisfies rjk = Ck, (md 
rji < ■ ■ ■ < rjn- Moreover, ryi — >■ as n ^ oo, L —>■ oo with n/L bounded above. 

We refer to the solution in Theorem 6.2 as the cut- ojf ground state (CGS). Our 
interest is the infinite length, finite density limit, and then the infinite density limit. 
For this purpose we recast (6.5) as a pair of integral equations. Assuming still that 
n = p, we let A = 27rn/L and define functions on the interval [0, A) by 

vio-vk, a(o = a M|^<^<2^. 

(We hope that context will prevent confusion between the use of r] and ( to denote 
functions and the use of 77 as C as pure variables.) 

Equations (6.2) with n = p are equivalent to integral equations 

-1 /•2'iTn/L 

viO = m + - / tan-i(r7(0 + CiO) d^, (6-10) 

Jo 

-I i>2nn/L 

It can be shown that for every set of data /, J corresponding to sets of integers 
indexed as in (6.6), there is a solution (77, C) with positive jumps at each of the points 
of discontinuity 2'nk/ L. We define approximate inverse functions accordingly: 

i(n) = : ^{i) > ry}; ^(C) - mf{^ : C(0 > C}- (6.11) 

Remark. These functions have jumps of size 2tt/L, so their difference quotients 
with respect to rj and to C, respectively count the number of particles associated to 
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With the Bethe equations in the form (6.10), it is clear how to pass to the infinite 
length, finite density limit. To pass also to the infinite density limit we observe that 
(6.10) is equivalent to 

77(0 - ^(0) = /(O - ^(0) (6.12) 
+ - / [tan-i(r/(e) + ((CO) - tan-i(r/(0) + ((C))] ^C; 

C(0 - C(o) = J(0 - J(o) 

-1 /•27rn/I/ 

+ - / [tan-i(r/(e') + ((0) - tan-i(r/(^') + C(0))] c/C'- 

The data (6.8) give rise to the functions = J{^) = ^ — ca in the L ^ oo limit, 
A finite, and in this limit ?7(0) = ^(0) = 0. Thus in the form (6.12) we pass to the 
A = oo limit and find the equations for the ground state {rjo, (q — fjo)'- 

^o(e) = e + - / [tan-i(r7o(e) + Voi^) - tan-\r,oim d^- (6-13) 

We obtain existence and uniqueness of the solution to (6.13) by converting it to 
a linear integral equation for the derivative of the inverse function: 

Diffierentiate (6.9), multiply by po, and change the variable of integration to obtain: 

Remarks. 1. In view of the remark after (6.9), po may be thought of as the particle 
density function for the ground state. 

2. Note that the kernel of the integral operator K satisfies 

1 r dr,' 1 

Therefore T has norm 1/2 as an operator in L°°(1R+) and (6.16) has a unique 
solution, obtainable by successive approximations. The solution is positive and 
bounded. 

In the next section we need both the analytic continuation po,t of p to the negative 
axis and a dual function pg^^ that is the analytic continuation of the solution to the 
same equation but with a change of sign: 

Poiv) = ^ + lf if('+%2 = 1 + KphivU > 0. (6.17) 
In terms of these functions we define functions and 

Uv)- f P^Mds-, Coiv)- r PoA')ds, fjen. (6.18) 
Jo Jo 
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Theorem 6.3. The extension po,t of the particle density function po and the dual 
function Pq ^ satisfy the estimates 

< poM < 1 < Pit < 2, e iR; (6.19) 
Po,t(^) = l + O(l^), plt{v) = 2 + 0{l/r]) as rj ^ oo; 
^VPoAv) = -'^ + 0{log\r]\/\r]\), pQt = l + 0{l/r]) as rj ^ -oo. 

The functions and have asymptotics 

Co{v) = V + 0{l), Co{v) = '^V + 0{l) as ry ^ +oo; (6.20) 
7r6(?7) = -1 + C>(log|?7|/|?7|), Q{v)^V + 0{l) as rj ^ -oo. 

In particular the function is a diffeomorphism of the line. Its inverse function 
rjo satisfies the integral equation (6.13) on 

Proof. It follows from (6.16) that the solution extends analytically. The kernel 
is positive, so < po < 1- The asymptotics of po and of are easily verified 
from (6.17). The argument for Pq ^ and are similar. Finally, (6.13) results from 
substituting (6.16) in the definition of and using podr]' = d^' to change the 
variable of integration. □ 

We conclude this section with the integral expressions of various quantities as- 
sociated with a state described by a pair of functions ry, Q. For a finite number of 
particles the momentum P and energy E are 

n n J /.27rn/L j /■27rn/L 

The particle number n = p itself can be expressed in terms of the inverse functions: 

n= ^supe(r/) = ^supe(C)- (6-22) 
Therefore in the infinite length limit the particle densities per unit length are 

^sup^v) = P{v)dm ^supe(C) = (TiOdC. (6.23) 



7. FINITE PARTICLE EXCITED STATES 

In this section we construct particle-hole excitations of the ground state. These 
have finite positive energy with respect to the ground state in the infinite density 
limit, confirming the stability of the ground state. 

We begin with a single particle/hole excited state for an r] particle, with momen- 
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state with a high particle density, so that there are large momentum particles. De- 
note the functions for this CGS by 77, C = V- We change the CGS data by choosing 
k such that rjk is close to rjh and changing Ik to Ik — s, where the integer s is chosen 
so that Ik — s is close to L^o{r]p) /2n. In the integral equation picture (6.10) the 
new data is given by functions I, J — I. Here I — I except on an interval of length 
27r/L located near — CoiVh), on which it differs from / by approximately — ^h, 
= Co(^p)- This change in data changes the rj solution by a bounded amount on 
that interval and otherwise changes both r]{^) and C(C) by 0{1/L). Therefore we 
write the solution {fj, () of the perturbed equation as 

m = v{o+ei{o/L+f{a c(o = c(e)+e2(o/^ (7.1) 

where ei vanishes on the critical interval and / is supported on the critical interval. 
With a similar notation for energy and momentum, the excitation momentum and 
energy are 

SP = P-P=1. j\e,iO - 62(0] + /; (7.2) 

SE = E-E=^j^ [c^eiiO- 0262(0] d^ + cif. 

Here we abuse notation and write / also for the value of fj—rj on the critical interval. 
With large L the integral equations (6.10) give 

f^^p-^h + - f [te.n-\rjh + f + vi^)) - ts^n'^rj^ + rj{e))] d^'; 

-^0 (7.3b) 

62(0 ~ 2[tan-i(r/,, + / + 7/(0) - tB^n'^rih + r/(0)] 

, 1 [ei(r)+e2(0](i0 

Here (7.3a) is valid outside the critical interval. 

In the infinite length, infinite density limit (6.13) and (7.3b) show that / = 
Vh — Vpj V ^ Vh on the critical interval, as desired. In this limit the remaining 
equations (7.3) become 

. . _l r Mv) + e2{v')]po{v')drj' _ 



62(^7) = ^ / 

Jo 



[ei{ri') + e2{v)] Poiv') dv' 



i + iv + V'? 

+ 2tan~-^(?7p + 77) - 2 i&n~^ {r]h -\- rj). 

We write 6^(77) for ej{^o{ri))po{r]): The integral equation (6.18) for po allows (7.4) 
to be rewritten: 



eiiv) = K{e2){r]); 



(7.5) 
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In the limit the excitation momentum and energy are 



SP 



SE= — 
2tt 



—a 



1 [°° 

^ [eiiv) - e2{v)] dr] + {rip - r]h); (7.6) 
1 

— j [-ciei(?7) + 0262(77)] dr] - ci{r]p - rjh) 
[ei (rj) + 62 iv)] drj + rjp-rjh^ 

[ei(??) -e2(??)](i?7-?7p + 77/i 



1 

2^ 



+ 

(The last part of (7.6) follows from the second part, since 2a = ci — C2 and 2co = 

Cl + C2.) 

The system (7.5) and the equations (7.6) thus characterize this single parti- 
cle/hole excitation of the ground state. Again the fact that K has norm 1/2 implies 
that (7.6) has, for each choice of rjh and rjp, a unique solution. Positivity of the 
kernel implies that the solution (ei, 62) is negative when rjp — rjh is negative, so the 
excitation energy is positive. We can be more precise. 

Theorem 7.1. The excited state momentum and the excitation energy can be ex- 
pressed in terms of the (extended) functions ^0, of Theorem 6.3. They have the 
form 

5P = Pp-Ph = P{vp) - P{vh) = Uvp) - UVh); (7.7) 
5E = Ep-Eh = E{r]p) - E{r]h) 

These functions have the asymptotics 

Ep - 2a6^l-^''l as Pp -00, (7.8) 
Eh ~ -ciPh as Ph +00. 

Proof. The inhomogeneous term for the pair {rjp, rjh) is the sum of the inhomoge- 
neous term for the pair (rjp, 0) and the inhomogeneous term for the pair {0,rjh) and 
is an odd function of its two arguments. Therefore it suffices to consider the case 
Vp ~ Vi Vh = 0. We make the fj dependence explicit by writing ei{rj,fj), e2{rj,fj). 
Set 

A^^|^le:+e.]; /, . ^|?[e, - e,]. (7.9) 

(Note that the Cj vanish at 77 = 0, so the sum and difi'erence can be recovered 
by integrating the fj from fj = 0.) Adding and subtracting first derivatives with 
respect to fj in (7.5), we find that as functions of rj the fj satisfy 

f, = Kf, + k, f^ = -Kf2 + k, kirj,fj) = --—^—-^. (7.10) 

TT 1 + (?7 + 77)^ 

Note that k is the kernel of the operator K. Let Fj denote the integral operator 
that has kernel fj and the identity operator. Then (7.10) is equivalent to the 
operator equations 
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The integral fjiv-iV) dr] is the same as the function Fj{l){fj), where 1 denotes 
the function that is identically 1. It follows from (7.9) and (7.10) that 



1 

2^ 



1 

2^ 



ei(?7, r])dr] + r] 



(/i-/2) + 2 



(7.12) 



= ^£[(i^i(i) + i) + (1-^^2(1))]; 



62(77,77) dry 



1 rv r /-oo 

2/ [I 



+ /2) 



(i?7 



= I J\f^{1) + F2(l)] dr7 = ^ + 1) - (1 - ^2(1))] cir7 

In view of (7.11), -Fi(l) + 1 is the solution of w = Ku + 1, while 1 — -^2(1) is the 
solution oi V = —Kv + 1. By (6.17) and (6.16) these functions are pg.t po, 
respectively. It follows that 

(7.13) 



1 f 

— J ei{v,v)dr] + fj=^[^*{f]) + ^o{v)], 
27r Jo ^) = ^ [^0 iv) - Co {fj)] . 



(7.14) 



The momentum term P{fj) in (7.6) is the difference of the two expressions (7.13) 
and (7.14). Similarly, the excitation energy term £^(77) in (7.6) is is the linear 
combination of (7.13) and (7.14) with coefficients — ci and C2 respectively. This 
proves (7.7). The asymptotics (7.8) follow from (7.7) and the asymptotics (6.20) of 
^0 and Q. □ 

Exactly similar considerations apply to a single particle/hole ( excitation, with 
the roles of 61 and 62 interchanged. In particular, the indices 1 and 2 in (7.9) should 
be interchanged throughout. 

Compound excitations may be treated in exactly the same way. Consider the 
case of particle/hole excitations in 77 described by r]h,i > 0, ?7p,i < 0, i = 1, . . . ,r 
and in ( described by (h,j > 0, (pj < 0, j — 1, . . . , s. Then in the infinite density 
limit the Bethe equations become 

s 

61 (r/) = K{e2)iv) + 2 ^[tan-i(r/ + (pj) - ta,n-\r] + (hj)], 

j=i (7.15) 

r 

62(77) = i^(ei)(?7) + 2^[tan"^(?7 + ?7p,i) - tan"^(?7 + 77^,^)]. 

This state is a superposition of corresponding single particle/hole excited states. 
Therefore the associated momentum and excitation energy are given by 



(7.16) 



sE = -a( J2ro{vp,i) - eo{vh,^)] + E[^o(Cp,i) - CoiChj)] 



i=l 



+ Co ( Yl^^o^^p,i) - ^oivh,i)] - Yl^^o^'^p,j) - Co(C/ij)] 
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8. FINITE DENSITY PERTURBATIONS OF THE GROUND STATE 

In this section we study large deviations from the ground state. Again we mo- 
tivate the procedure by looking briefly at the finite length, finite particle case in 
the form of the integral equations (6.10). Denote by Iq = Jq the data that corre- 
sponds to the CGS data (6.8). The CGS data uses all the integers from an interval 
with cndpoints ±0{n/2) and we are temporarily applying a low energy cut-off that 
corresponds to the upper end of this interval, so the data that corresponds to an 
admissible perturbation of the CGS satisfies I, J < Iq. The solution can be found 
by iteration from the integral equations (6.10). The monotonicity properties of the 
mapping imply that the solution is dominated by the CGS (r/OiCo)- 

v{0<vo{0, C(0<Co(0- (8.1) 

Moreover if the perturbation is to have finite relative energy density in the limit, 
we expect that /, J ~ /q for large ^. Note also that the jumps in the data must be 
at least as great as the jumps 2tt/L in Iq. We pass to the infinite length, infinite 
(total) density limit and look for a solution of an analogue of (6.12): 



1 r°° 

viO - m) + - / [tan-^(r7(e) + CiO) - te.n-\r,o{0) + Co(e'))] ^C; 
m - AO + - / [tan-i(r7(r) + C(0) " t^^-'iVoiO + Co(0)] d^- 

71" Jo 



(8.2) 



According to the preceding discussion the natural assumptions on the data are 

^>i, ^>i; m)^Jio-^^, e»o. (8.3) 

It can be shown that finite particle solutions whose data is such that the formal 
limits are (8.2) converge to solutions of (8.2). However it is more convenient for 
what follows to write the equations for the particle density functions p — d^/drj, 

,(,) . nn) = 1 - i £ = 1 - (K.)(r,y, (8.4) 

^(C)-'(C) = l-i/_:,^ = l-(.p)(C). 



where 



We have removed the upper bound on particle energies, so the integrations in (8.4) 
need to be taken over the whole line. Now K has norm 1 as an operator in L°°(]R) 
so solvability of (8.4) becomes an issue. We rewrite (8.4) in the form 

p = e{l-Ka), a = -^{l-Kp) (8.4') 



and adopt the point of view that the functions and \1/ are the data. For the 
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s > 0. More generally we can express the conditions (8.3) and the assumption that 
our states are perturbations of the ground state by requiring 

< e, * < 1; lim (e - Bo) = lim (* - *o) = 0, (8.5) 

|??|-<-oo ICHoo 

while (8.1) gives an implicit constraint on the data (G, \1/) in terms of the solution. 
In fact (8.1) implies the reverse relation between the inverse functions, leading to 

(P-Po)>0, f ((T-(To)>0, aUr7,(7. (8.6) 

Moreover the cut-off solution (ry, () and the corresponding CGS (?70) Co) at any given 
step were defined over the same ^ interval, so the relative particle density should 
be zero: 

/oo pOO 
(p-po) = 0= / (a-ao). (8.7) 
-oo J —oo 

We show below that (8.7) is necessary for the relative energy density to be finite. 

Theorem 8.1. Under assumption (8.5) the equations (8.4) have a unique solution. 
For the solution we have 

0<p<l, 0<(J<1; (8.8) 
lim p{r]) = lim a{C) = 1; lim p{r]) = lim a{C) = 0. 

TJ^OO >00 TJ— > — 00 >— 00 

Proof. We use and \1/ also to denote the operations of multiplication by the 
corresponding functions. Then the equations are 

p + OKa = 0, a+ ^Kp = 

We apply the operator Id — to the first equation and Id — QK to the second 
and reorganize to obtain 

{Id - QK^K)p = - 0K*, {Id - ^KQK)a = * - ^KQ. (8.9) 

It is enough to show that QK^ and ^KQ have norm < 1. The assumptions on 

and \E' imply easily that the integral of the kernel of either operator against 
either variable is bounded away from 1, so the operators have norm < 1 in any L^, 

1 < p < 00. The asymptotics follow easily from the assumptions about and 
together with the trivial estimates 

O<p<0; 0<a<^. □ 

One must make stronger assumptions to ensure finiteness of the integrals (8.6) 
and (8.7) and the relative energy density 

E = -^J^{p- po)vdv +^J^(^- ^o)CdC. (8.10) 
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Theorem 8.2. The integrals (8.6), (8.7) are finite if (6 - Gq) and - ^'o) are 
integrable. 

If vi^o — 0) cind C(^o — ^) ore integrable, then E is finite if and only if (8.7) 
holds. If so, then E > 0. 

Proof Note that 

p - po = (e - eo)(l - Kao) - QK{a - ao), (8.11) 

with a similar equation for a — ao. Analyzing this system as above, we see that 
integrability of p — po cind a — ao follows from integrability of the data in the 
analogue of (8.9), which follows in turn from integrability of — ©o and — ^^o- 
Note that integrability of r]{po — p) at —oo follows immediately from (8.11) and 
integrability of ?7(©o ~ ©)• At -\-oo we observe from (8.11) that the necessary and 
sufficient condition is the integrability of r}K{a — ao)- Assume that C^(^ ~ ^o) 
is integrable at — oo. Then (^"^{a — ao) is also integrable at — oo and Lebesgue's 
Dominated Convergence Theorem implies that 

/oo 
{a-ao)dC,. (8.13) 
-oo 

This proves the necessity of the conditions (8.7). 

Conversely one can deduce the integrability of r]K{a — ao) and C^iP ~ Po) from 
conditions (8.7) together with integrability of ri{po — p) and ({ao — cr). To see this, 
note that the mean-value zero condition allows us to replace the kernel of K with 
the kernel 

1 1 



%C) = - 

TT 



1 + (^ + C)^ 1 + '. 



(8.14) 



The corresponding operator takes L^(R, (1 + 1^1) d^) to L^(R, (1 + \r]\)dr]), since 
we can multiply k hy {1 + \ri\) /{I + {z])., take the absolute value, and estimate the 
integral with respect to rj. 

Finally, to show positivity of the relative energy density we let gi and g2 denote 
the integrals in (8.6), so that 

^ = p, ^ = a; gj>0, gj{±oo) = 0. (8.15) 

Integrating by parts in (8.10) we find that 



9. THERMODYNAMICS 



We consider now the thermodynamics of the finite density perturbations of the 
ground state that were constructed in the preceding section. Wc follow, with ap- 
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finite L, a state with a very large number of particles, and data I, J as in (8.1). 
Define 

h{s) = s-- r[tan-i(s + C(e'))-tan-i(^(0) + C(e'))]^e', ^ > ^(0); 

(9.1) 

m=t-- r[tan-i(r7(O+^)-tan-i(r7(O+C(0))](ie', t > C(0). 

These are strictly increasing functions with the properties that at the (discrete) 
values taken by the data the values are 

%(0) = no, Kcm-Jio- (9.2) 

The number m oi rj particles from our state with momenta in a small interval 
(?7, rj + drj) is proportional to the length of the corresponding ^ interval, i.e. if 
V = ^(0: V + dv ^ Vii + then 

m ^d^. (9.3) 
The function / has jumps that are positive integer multiples of L/27r, so 

m<n=^(/(e + rfO-/(0)- 

If data / at the m jump sites in the interval (^, ^ + d^) is changed in such a way 
that / remains strictly increasing and the values of LI/ 2% are integers, the effect 
on the total energy of the corresponding state will be negligible: 0{{dr])'^). The 
possible sites are indexed by the n integers between L/(^)/27r and L(/^)/27r. Thus 
the contribution to entropy from the interval is log ( ) . 

In the infinite length and infinite density limits the density functions for the 
numbers n and m above have densities with respect to the momentum rj given by 
p = d^/dr] and pt = dh/drj respectively. As before p is the particle density. We take 
Ph = pt — P be the hole density and pt to be the (formal) total density. Similar 
considerations apply to ( particles and give particle density a and hole density 
o'h = CTt — CT- In the notation of §8 the relations among these functions are 

p = Qpt, Ph = {l- Q)pt; o- = *(7t, ah = {l- *)(7t; (9.2) 

pt^l- Ka ^1- K<^au at = l- Kp = l- KQpt. 

Returning to the interval (^, ^ + d^), let (?], i] + drj) denote its image under h~^. 
According to the preceding argument and Stirling's formula, the entropy density 
over {r] + drj) is approximately 

^[pt^ogpt- plogp- phlogph] (9.3) 
= :^[p^og{pt/p) + Phlog{pt/ph)] 
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A similar formula holds for particles of the second type. Thus in the infinite length, 
infinite density limit the entropy per unit length is 

1 r°° 

S =— [ptlogpt- plogp- phlogph]dri (9.4) 

1 f°° 

27r J (Tt-aloga- an log an] dC. 

The entropy density for the ground state is zero. 

At a given finite temperature T the free energy density of a perturbation (p, a) 
of the ground state relative to that of the ground state is 

Ft{p,(t) = E-TS (9.5) 

where the relative energy density E is given by (8.9). In order for E to be finite, 
we showed that the relative particle densities must vanish: 

/oo poo 
{p-po)dr] = 0= / {a-ao)dC. (9.6) 
-oo J — oo 

Thus the equilibrium solution should minimize F subject to the constraints (9.6). 

We consider briefly the formal argument for an equilibrium. Equations (9.2) 
relate variations in pt, ph, crt, and ah to variations in the particle densities p, 0: 

Pt = -K&, pf^ = -p- Ka; (9.7) 
at = -Kp, ah^ -a - Kp. 

The constraints (9.6) require 

/oo j-OO 
pdr] = 0= / adC- (9.8) 
-oo J — oo 

Taking into account (9.7) and the symmetry of the operator K, we can express the 
variation in the free energy as 

/oo 
p[-cir] + T \og{p/ph) + TK \og{at/ah)\ drj (9.9) 
-oo 

/oo 
a [c2C + Tlog{a/ah) + TK\og{pt/ Ph)] dC,. 
-oo 

Following Yang and Yang^, one can introduce functions e^: 

e-^i/^ = p/Ph, e-'-l^ = a/an. (9.10) 
and derive Euler-Lagrange equations in the form 

ei = AiT - ciT] + TK log(l + e'^^/^), (9.11) 

, A rr \ „ ^ \ rrT^^ — I „ — ei/T\ 
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where the Aj are the Langrange multiphers that correspond to (9.8). 

As noted in the introduction, there are two difficulties with this argument. First, 
equations (9.11) have no solution: taking derivatives one sees that the ej should be 
decreasing functions, so the functions log(l + e~'^^/^) are increasing and bounded 
above, but they are taken by K to decreasing functions and thus one needs ei < 
— ci?7 + 0(1) at +00. Therefore log(l + e~'^i/^) is nonnegative and grows like ry at 
+00, so K\og{l + e~^'^/'^) cannot be defined. Second, it is not clear that the free 
energy is bounded below. 

The first difficulty is easily dealt with. We may assume that r and s are test func- 
tions, and (9.9) simply means that they are derivatives of test functions. Therefore 
the condition F — simply means that the derivatives of the two bracketed terms 
vanish. Now K anti-commutes with derivation, so the correct form of (9.11) is 



dei 
dr] 



-ci + K ^' 



de2 



Combining this with 



G = 



-6l/T 



1 + 



* = 



-62 /T 



1 + e-e2/r 



(9.13) 



yields a well-behaved system of equations for the thermal equilibrium particle dis- 
tribution (p, a) = {p^^\a^^^). 

To deal with the second difficulty we prove the following. 

Theorem 9.1. The free energy functional Ft is convex and bounded below by 
—CT^, where C is a constant that depends only on the velocities Cj. 

Proof. From (9.9), the second derivative of the free energy for the variation (p, a) 
is 



oo 

Tl p 

-OO 

+ 



= T 



+ T 



Ph 

"OO 



drj 



,/_oo [O- (Th \pt Ph 

+ {K&) 



P Ph 



(Th 



(^hJ \(Tt CTh 



(9.14) 



(- 


_ ph_\ 


\pt 


Ph) . 



dr] 



dz. 



We use (9.7) and reorganize to find that the second derivative is 



\VPPh \fPtPh 



Kp\ dC (9.15) 



To prove boundedness we introduce the function 



It 



,^ 1, 
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Then G(6) = GdGo — ©1) and the entropy density can be written 

PtG{e)dr]+ / (7tG{^)dC. (9.16) 

-oo J —oo 

For the first integral we partition the hne into three sets 



h = {v 
h = {n 



G(e)<|?7(eo- 6)1/5}, (9.17) 
|0o-0| > 1/2, G(e) > 177(00-0)1/5}, 
|0o - 0| < 1/2, G(0) > |r7(0o - 0)|/5}, 



where the constant 5 is to be chosen. Note that < G{t) < log 2 so r/ G I2 implies 
|?7| < 25 log 2. Also 

-logt> -2tlogt>G{t), 0<t<l/2 

so ?7 e /a implies |0o — *| < 6xp(— 1?7|/5). Now —tlogt is increasing on the interval 
[0, 1/e], so we conclude that r] E I3 implies that either \r]\/5 < 1 or 



^(0) = ^(100 - 01) < -2100 - 01 log 100 - 01 < 2Me-l'?l/^. 



Since 1 < 2 log 2 we can combine these estimates and obtain 

/oo 1 ^cxj POO I I 

PtG{e)dr]<- pi(0o-0)?7rf?? + 451og2 + 2 / ^e-l'^l/^dr] 
-00 <^ J-00 J -00 <^ « 

/oo 
Pt{eQ-e)r]dr] + mog2 + U. 
-00 

We claim that 

/oo 1*00 1*00 

Pt{Qo-Q)vdv< {po-p)vdv+ {ao-a)CdC. (9.19) 
-00 J —00 J —00 

In fact (8.11) can be rewritten 

pt(0o - 0) = p - po + 0o^(p - Po). (9.20) 

Using ( , ) to denote the standard inner product, we want to prove that 

/oo 
0oK(a - ao)ri dr] = {K{a - ao) , Qov) (9-21) 
-00 

< (o-o - cr , C) = -{Dg2 , z) = {g2 , 1), 

where D denotes differentiation and (72 is the function in (8.15). We use the fact 
that derivation anti-commutes with K and integrate by parts to obtain 

{K{a - ao) , 00^7) = {KDg^ , 0or?) = -{DKg^ , 00^7) 



(9.18) 



30 



Now < KOo < 1 and g2 > 0, so this proves (9.19). Using (9.18), (9.19), and the 
analogous results for cr, we obtain 

Ft{p, (t) < (ci - p (po - P)vdv + 45(1 + log2) (9.22) 
+ f -C2-^j J (ao-a)CrfC + 45(l + log2). 

We take 5 = 2T/(ci - C2) and obtain C = 8(1 + log2)/(ci - C2). □ 
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